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Abstract 

Motivated by the forthcoming data from the E950 experiment at BNL for small an- 
gle polarized proton-carbon scattering we revisit the problem of Coulomb phase. The 
approximation usually used, with the momentum transfer squared (q 2 ) small com- 
pared to the inverse elastic slope, is not justified within the kinematics of the E950 
experiment. We go beyond this approximation and derive a new rather simple expres- 
sion which recovers the formula of Bethe at small q 2 , but is valid at any momentum 
transfer. 
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The forthcoming precise data for small angle polarized proton-carbon scattering mea- 
sured in the E950 experiment at BNL raise again the problem of Coulomb phase QXJ] . The 
experiment covers the range of momentum transfer up to q 2 = 0.05 GeV 2 where the product 
Bq 2 = 3 (B « 60 GeV 2 is the elastic slope). At the same time the early calculations of the 
Coulomb phase neglected terms of the order of Bq 2 [l], [|. The next-to- leading corrections 
of the order of q 2 ln(g 2 ) were calculated by Cahn ||. However, next-to-next-to-leading order 
terms might be also important at higher q 2 . 

The asymmetry of polarized proton-nucleus scattering arising from Coulomb-nuclear 
interference [Q, the main goal of experiment E950, can be essentially affected by the Coulomb 
phase, especially if the hadronic spin-flip amplitude has a real part []5J. This situation 
motivates us to attack once again the problem of the Coulomb phase to reveal the corrections 
of the order of q 2 and higher. This has recently been done numerically for pp scattering in 

We demonstrate that one can arrive at a rather simple analytical expression Eq. (|34]) 
for the Coulomb phase which is valid at any q 2 , provided that the electromagnetic formfac- 
tor and nuclear amplitude have Gaussian dependence on the momentum transfer q. This 
assumption is well justified for light and medium heavy nuclei. For a proton target the 
Gaussian formfactor falls off too steeply at q 2 > A 2 , where A 2 = 0.71GeV 2 is the pa- 
rameter of the dipole parameterization. However, the Coulomb phase matters only in the 
Coulomb-nuclear interference region q 2 <C 0.1 GeV 2 where the Gaussian parameterization 
is precise. 

If the two slopes of q 2 dependence for electromagnetic and hadronic amplitude coincide 
the expression for the Coulomb phase becomes especially simple as is given by Eq. (|35|) . 
The elastic scattering amplitude can be represented as, 

f(q) = Mq) + f N c(q) , (1) 

where the first and second terms correspond to the first (a) and sum of the second and third 
(b,c) graphs depicted in Fig. Ill respectively. Of course such a classification is conventional 
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Figure 1: Three types of interaction: pure Coulomb (a), nuclear (b) and 
nuclear-Coulomb (c). 



and we make it for convenience. 
The amplitudes in Eq. ([!]) read, 

fc(q) = ^Jd 2 be^ b (l-e^) , (2) 
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fNc{q) = ^jd 2 be* b e i ^ lN (b) , 

where 

lN(b) = ^- [d 2 qe-^ b f N (q) . 
The amplitudes Eqs. (P - (|) are normalized according to the conditions, 

da 
dq 2 

Vtot 
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/(g) 

Anlmf N (q = 0) . 
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In Eqs. (0) - fl5|) fc(q) and /jv(g) are the net Coulomb (long-range) and nuclear (proton- 
nucleus) amplitudes (Figs. |l|a and b respectively), and fNc{q) includes the effects of strong 
interaction (Fig. |l|b) and Coulomb-nuclear interference (Fig. |]c). 

The Coulomb phase Xcip) m Eqs. (|I|)-(^) is a Fourier transform of the Coulomb part of 
the amplitude calculated in the Born approximation, 
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where F^f'(q 2 ) are the electromagnetic formfactors of the colliding particles (nuclei), and 
Zip are their charges. In order to keep the integrals finite we give to the photon a small 
mass A which will disappear from the final expressions in the limit A — ► 0. 

On the contrary to the Born approximation, the full Coulomb amplitude including all 
the multi-photon exchanges is complex, i.e. has a phase <&c{q), 



fc(q) = -sign(Zi Z 2 ) f c (q) 



where 



fc(q)\ = \f { c B) (q)\ [l + 0((aZiZ 



(9) 
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For the sake of simplicity we restrict our consideration to the case of a | Z1Z2 | <C 1 
(appropriate for pC) and will neglect the higher order corrections. In this approximation, 
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2 jd 2 b X c{b)e^ b ' 

' ' d 2 qi d 2 q 2 f { c B \qi) ff\q 2 ) 5{q - qi - q 2 ) . 
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2*7c (?) J 

With the same accuracy the amplitude fNc(q) in Eq. @ can be represented in the form, 

fso(q) = f N (q) e^ cM + 0{{aZiZ 2 ) 2 ) , (12) 

where ^ 

fNc(q) $Nc(q) = 7T d 2 qi d 2 q 2 5(q- qi - q 2 ) . (13) 

ZTT J 



Our results for $c, &nc and A<3> = $c — &nc look rather simple if the g-dependence of 
the formfactors in Eq. (|3|) is Gaussian, 



-a q 



f N (q) = f N (0)e 



-bq 2 



«=g((r 2 )i + (r 2 ) 2 ) 
b = B/2 , 



where B is the slope of the hadronic differential cross section. 
In this case the phase Eq. (|13|) takes the form, 



<&Nc{q) 



aZiZ 2 f d 2 q\ 
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q{ + A ; 



exp 



-(a + b) qf + 2bq 1 -q 



This integration can be performed analytically replacing 

1 



q\ + 

Then the phase Eq. ( jlBf) takes the form, 

a.Z\Zi f du 



J dt exp -t (q 2 + A 2 
o 



®Nc(q) 
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where 



b 2 q 2 
a + b 

A 2 (a + b) . 
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Z = 
V = 

Further, expanding the exponential we get 

<5>Nc{q) = 
where 

E k +i{v) 
Ei (v) \ v ^ 

Here 7 = 0.5772 is the Euler constant. 
Thus at v << 1 we arrive at, 



aZ 1 Z 2 ™z k 

2^ 7T E k+ i(v) e , 
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k > 
-7 — ln(i>) . 



®Nc(q) = olZ x Z 2 [27 + huW) -Ei(z)\ +0(v) 
where Ei(z) = — Ej(— z) is the integral exponential function. 



(14) 
(15) 
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(17) 



(1* 
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(22) 
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Coming back to the phase ( |TTD we represent the integral in the numerator in the r.h.s. 



as, 



U'D = ^ / Al A2 fc\'h) /<?'(*) S(q- ?i - ? 2 ) 



2A 



q + A 



f(B) 



g - A 



(23) 



Since the product 



S 



q + A\ s ( q- A 



exp 



A 



(24) 



is independent of the angle between the vectors g and A (d 2 A = dA 2 d<f>/2), one can 
perform integration over in (^3), 

OO 2-7T 



ie(?) 



2 (aZxZa) 
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exp 
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x [(g 2 + 2gA cos <\> + A 2 + 4A 2 ) (g 2 - 2gA cos + A 2 + 4A 2 ) 
Then, it can be represented as, 

OO 

J c (g) = 2 (aZxZ 2 ) 2 / dA 2 exp 



-i(g 2 + A 2 ) 



^(A 2 ,g 2 ,A 2 ^ 



where 



^(A 2 ,g 2 ,A 2 ^ 
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q 2 + 4A 2 



g 2 + A 2 + 4A 2 ' ^ 

In the limit A — * function \1/(A 2 , g 2 , A 2 ) can be expanded dependent on the relation 
between A 2 and g 2 , 
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Performing integration in (|26"[) in parts and taking into account relation (|29|) we arrive 



at, 



h{q) = f { c B \q)^ c {q) = -aZ 1 Z 2 f { c B \q) 

'2 + t 



/ 2\ 

In I |jj +w y die- 1 " 1 In 



2 + t 



w / die - "" In 



+ 0|- 



A 2 ^ 
4 2 



(30) 



5 



where 



w 



a q 



(31) 



The integrals contained in this expression can be represented in terms of the integral expo- 
nential function, 



o 

w J dt e~ wt In 
-l 

w I dte- wt In 



2 + t 
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2 + t 
t 



Ei(w) - 7 - \n(2w) + e 2w fEi(io) - Ei(2iw) 
7 + ln(2w) + e 2w E 1 (2w) . 



(32) 



Eventually we arrive at the following expression for the phase <&c{q), 



<S>c{q) = aZ x Z 2 I 2 \n(2w) - In ^ +2 7 - Ei(w) + e 2w [2 Ei(2w) - E x (w 



(33) 



Apparently, the common phase factor for the terms in the elastic amplitude (|l|) is unob- 
served, only the difference A$ = $c — $>nc, the so called Coulomb phase, matters, 



A$(g) = aZ x Z 2 \ In 



F +Ei(z) 



Ei(w) + e 2w 2Ei(2w) -Ei(tu) 



(34) 



Note that all the divergences and the photon mass A have cancelled. 

If the slopes of the Coulomb and hadronic formfactors coincide, a 
the expression fl3"lp takes an especially simple form, 



A$(g) = aZ x Z 2 e 2w 2E 1 (2w)-E 1 (w 



b (z = aq 2 /2 = w), 



(35) 



This is a good approximation for nuclei where the slopes of the electromagnetic and hadronic 
formfactors are mainly determined by the nuclear radius. For a proton target the electromag- 
netic slope a = 4/A 2 = 5.6 GeV~ 2 . The slope of the hadronic amplitude b is approximately 
equal to a at medium-high energies. However, it increases with energy and at very high 
energies, in particular at the Tevatron, b substantially exceeds a and one should use the 
exact expression (|34T). 

To see the difference of our results from the previous calculations we show A$(g 2 ) by the 
solid curve in Fig. |2|. We performed calculations for proton-carbon elastic scattering using 
(|35|) with a = 30GeV~ 2 . Dotted and dashed curves correspond to the formulas of Bethe 
Jl| and Cahn || (with B = 8/ A 2 = 2a) respectively. At small q 2 dominated by the leading 
~ ln(g 2 ) term the results of Cahn and ours coincide and deviate from the Bethe's due to the 
next-to-leading corrections ~ g 2 ln(g 2 ). At higher q 2 the next-to-next-to-leading corrections 
~ q 2 become important and all three curves are quite different. This is not surprising since 
it is not legitimate to use the results of |TJ, |3[] at so large q 2 . 

Note that the eikonal approximation we use is subject to inelastic corrections. In addition 
to the ordinary inelastic corrections to the hadronic amplitude which we assume to be 
already included, excitation of inelastic states N* between the electromagnetic and hadronic 
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Figure 2: The Coulomb phase as function of q 2 . The dotted, dashed 
and solid curves corresponds to the formulas of Bethe jj^, Cahn ^j, and 
present calculations, respectively. 



vertices in Fig. [I|c can affect our results. Nevertheless, it turns out that the correction to 
the Coulomb phase is quite small. Indeed, the N — > N* amplitude contains an extra factor 
q compared to the elastic iV — > N one. At small q —>■ this vertex squared cancels 1/q 2 
in the photon propagator and this correction looks similar to the ordinary inelastic one. 
We only have to replace factor A 2 jR\ in the hadronic inelastic correction to AZot/R 2 A . 
Then, aZ is the common factor for the Coulomb phase, and the relative correction becomes 
l/[(m%* —m%)R\]. This is a small correction, about 10~ 2 , to the constant and q 2 R\ terms 
in the phase. 

Concluding, we derived a new expression (|34j) for the Coulomb phase which is valid 
at any q 2 provided that the electromagnetic formfactor and the hadronic amplitude have 
Gaussian dependences on q. If their slopes are equal the Coulomb phase takes a very simple 
form fl35l) . 
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